On the zero temperature critical point in heavy fermions 
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I. INTRODUCTION 



Many of the physical properties of heavy fermions are 
determined by the proximity of these systems to a zero 
temperature critical point p]]. The existence of a zero 
temperature phase transition, from a phase with long 
range magnetic order to a non-magnetic state, allows for 
a formulation of a scaling theory which gives a very sat- 
isfactory description of these materials Qj. This scaling 
theory is most easily formulated within the context of the 
Kondo lattice model which emphasizes the importance 
of spin fluctuations neglecting charge fluctuations. This 
model has two relevant energy scales, namely the band- 
width of the large conduction band W, and the interac- 
tion J between the local f — electrons and the conduction 
electrons in the wide conduction band. At T = the rele- 
vant quantity is the ratio J/W and the zero temperature 
phase transition occurs at ( J/W) c which is different from 
zero for d > dr,, the lower critical dimension (probably 
dr, = 2). For finite temperatures the phase diagram for 
d > dr,, where a long range ordered magnetic phase exits 
at finite temperatures, is shown in Figure 1. 

This magnetic phase, generally of the antiferromag- 
netic type occurs for (J/W) < (J/W) c and below a crit- 
ical Neel line. 

For (J/W) > (J/W) c and finite temperatures we 
pointed out the existence of a crossover line, which we 
identified as the coherence line and marks the onset, with 
decreasing temperature, of a renormalized Fermi-liquid 
regime on this non-critical part of the phase diagram jlj . 
This physical interpretation of the crossover line allows 
to recognize immeaditely that a system just at the criti- 
cal point, i.e., with (J/W) — (J/W) c , should behave as a 
non- Fermi liquid since it does not cross the coherence line 
and consequently does not enter the Fermi liquid (FL) re- 
gion of the phase diagram In this paper we obtain 
the thermodynamic properties of this critical system. We 
also discuss and clarify the question of breakdown of hy- 
perscaling below the coherence line. 
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FIG. 1. The phase diagram of heavy fermions close to the 
zero temperature critical point (schematic). In the present 
case the shift exponent of the Neel line ip = 1 is different from 
the crossover exponent vz > 1 characterizing the coherence 
line T c . The crossover lines T c \ and T C 2 are both governed 
by the same shift exponent of the Neel line Below T c 
hyperscaling is violated and the thermodynamic properties 
are governed by mean-field exponents. 
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II. RENORMALIZATION GROUP EQUATIONS 



III. GENERAL SCALING RELATIONS 



The scaling results and the equations for the Neel and 
coherence lines can be obtained considering the expan- 
sion of renormalization group (RG) equations close to the 
zero temperature fixed point at (J/W) c . These expan- 
sions have a general character and can be written as, 



K n+1 =K C + b x ' v (K n - K c ) - (T/J) n 



(1) 



(2) 



The first equation describes the renormalization of the ra- 
tio K = (J/W) close to the critical point K c — (J/W) c 
under a change in length scale by a factor b. The expo- 
nent v is the correlation length exponent. Eq. 2 describes 
the renormalization of the temperature close to the zero 
temperature critical point. Since temperature is a pa- 
rameter its scaling is controlled by that of the coupling 
constant J which at the critical point scales according to 
the dynamical exponent z [0]. In Eq. 1, temperature en- 
ters linearly which corresponds to the lowest order term 
in an analytic expansion in powers of (T/J). The solu- 
tion of these recursion relations to a length scale I = b n 
is given by: 

Ki=K c + I 1 1" (K - K c - a(T/J)) + al z (T/J) (3) 
where a is a constant. Since I is arbitrary we iterate until 

l l ' v (K - K c - a{T/J)) = 1 
and this length scale defines the correlation length 



l = i = \K -K c -a(T/J)\ 



\3- join- 



where j = (J/W) and j c (T) = (J/W) c + a(T/J). At this 
length scale, substituting I for £ in Eq.3, we obtain 



K { = K c + 1 



a(T/J) 

\j-un 



The line T c /J = \j — j c (T)\^ z in the non-critical side of 
the phase diagram is a crossover line since K j and con- 
sequently all functions of this ratio have different asymp- 
totic behaviors for T >> T c and T «T C . It is governed 
by the exponent vz and is identified as the coherence line 
in Figure 1. On the other hand, at the line j — j c {T) or 
Tjm I 'J — ~\j — j c (T — 0)| the correlation length diverges 
and this line can be identified with the Neel line in the 
critical region, i.e., for (J/W) < (J/W) c . If we define a 
shift exponent if) through the relation Tj^/J oc \j — j c \ 
[^],[|, the present approach leads to ip = 1 independent 
of the crossover exponent vz. The critical line vanishes 
linearly close to the zero temperature critical point which 
is due to the fact that in Eq.l the lowest order analytic 
term in temperature is a linear contribution. 



The general recursion relations, Eqs.l and 2 give rise 
to the following scaling form for the singular part of the 
free energy density close to the zero temperature critical 
point § " 



fK\5\ 2 - a f [T/T c , 



(4) 



where S = j — j c (T = 0) and T c oc \S\ UZ . In this section 
we consider the scaling properties of physical quantities 
in the non-critical side of the phase diagram, i.e., for 
(J/W) > (J/W) c . We shall neglect the temperature 
dependence of j c (T) which is valid as long as vz > if) 
and T < T c . At T = , /o(0) = constant and the above 
expression gives the singular part of the ground state free 
energy density close to the transition (8 = defines the 
critical point). In fact at T = 0, Eq.4 defines the critical 
exponent a. The exponents a, v and z defined above 
are associated with the zero temperature unstable fixed 
point governing the quantum phase transition ffl . 

The existence of a FL regime for T « T c , as found 
experimentally, allows to take a Sommerfeld type of ex- 
pansion for the scaling function of the free energy density 
in this temperature range . We get 



f^\5\ 2 - a {a 1 +a 2 (^f + a 3 (^) i + 



(5) 



where a\, a 2 and 03 are constants. The specific heat for 
T << T c is then given by 



5T2 K 11 



T 



(G) 



which shows the characteristic Fermi-liquid linear tem- 
perature dependence with an effective thermal mass 
m T ocC/Toc \8\ 2 - a - 2vz . 

The scaling expression for the free energy can be easily 
generalized to take into account the presence of a field 
conjugate to the order parameter, in the present case a 
staggered field H and of a uniform magnetic field h [0, 
we get 



f^\S\ 



H 



\sr'\s\ 



/3+7 : 



\s\ q 



(7) 



where v and z are the correlation length and dynamic 
exponents as before, (3 and 7 are standard critical ex- 
ponents which describe the critical behavior of the stag- 
gered magnetization and susceptibility respectively. 4>h 
is a new exponent associated with the uniform magnetic 
field h. These exponents obey standard scaling relations 
as a + 2(3 + 7 = 2 but the hyperscaling relation is mod- 
ified due to the quantum character of the critical point 
and is now given by 



a = v(d + z) 



(8) 
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where d is the dimensionality of the system pL From 
Eq.7 we can deduce the scaling properties of the 
physical quantities of interest like the low tempera- 
ture (T << T c ) field-dependent uniform magnetization 



relation i.e. a 



i-(d+z) 



m{h) oc (df/dh) oc \6\' 



f m (h/h c ) with h c oc |<5| 



<Ph 



the uniform, low field, susceptibility, %o oc d 2 f/dh 2 
\5\ 2 - a - 2 * h fs{T/T c ), etc.. 

The scaling theory formulated above is completely gen- 
eral and any description of heavy fermions in terms of a 
zero temperature phase transition should fall within this 
approach. The crucial question is of course to deter- 
mine the critical exponents, i.e., the universality class of 
the zero temperature critical point. It is interesting to 
point out that for d = 1 < rft, (J/W) c = 0, but there 

is still a crossover line for (J/W) > which is given 

-i 

by ksTc oc We J / w , an expression similar to that of the 
Kondo line El . 



IV. GAUSSIAN THEORIES 

Eq. 6 for the specific heat in the Fermi liquid regime 
(T << T c ) can be written, using the modified hyperscal- 
ing relation, Eq.8, as 



m T oc 7 = C/T oc \8\ u{d ~ z) 



(9) 



For dimension d = 3 and in the case the dynamic ex- 
ponent z = 3, which is generally associated with ferro- 
magnetic fluctuations {| this gives essentially the result 
7 oc Ln|5| although the scaling theory cannot reproduce 
the logarithmic singularity. This logarithmic dependence 
has the same origin of the paramagnon mass enhance- 
ment of a nearly ferromagnetic system obtained origi- 
nally by Doniach and Engelsberg H] for the Hubbard 
model. In this case \S\ = \{U/W) - (U/W) c \ where U is 
the Coulomb repulsion, W the bandwidth and (U /W) c 
is the critical ratio Q| at which the T = ferromagnetic 
transition occurs. 

The case of antiferromagnetic fluctuations character- 
ized by z = 2 can be considered going to higher order 
in the Sommerfeld expansion of the free energy. Again 
using hyperscaling, i.e., Eq.8, we get 



7 = C/T oc a 2 \8\" (d - z) + a^\5\ u{d - Zz) T 2 



In three dimensions with z = 2 



this yields 



C/T oc a 2 \S\ 1/2 +a 3 



T 2 



3/2 



(10) 



(11) 



for T << T c where we took the Gaussian exponent 
v = 1/2. This result coincides with that of the Gaus- 
sian theory of Millis || for a nearly antiferromagnetic 
Fermi liquid and relies on the validity of the hyperscal- 
ing relation, Eq.8, below T c . In fact the exponent a in 
Gaussian theories of zero temperature magnetic phase 
transitions is determined by the modified hyperscaling 



2 -, where we used the Gaussian 
value for the correlation length exponent v = 1/2. Con- 
sequently in these theories hyperscaling is automatically 
satisfied. Notice that for d + z > 4 the Gaussian expo- 
nent a < implying that singularities in these theories 
are weaker than those obtained in a mean-field approach 
where the exponent a is fixed at a — for d + z > d c . 
Here d c = 4 is the upper critical dimension for these mag- 
netic transitions [QJ^] . However the mean field exponents 
a = and v = 1/2 violate hyperscaling above d c , i.e., do 
not satisfy Eq.8 for d + z > 4. 

Let us consider again the enhancement of the thermal 
mass for an itinerant three dimensional nearly ferromag- 
netic system on the light of the discussion above. In 
the general scaling approach we have, mr oc \5\ 2 a 2l/z , 
which is obtained without using the hyperscaling rela- 
tion. Using the mean-field exponents a = 0, v = 1/2 
and z = 3, which violate Eq.8 for d — 3, we get 
m T oc \{U/W) - {U/W)^ 1 = \8\~~ 1 . 

Since the mean- field susceptibility exponent j = 1, we 
find Xo oc \S\ where \o is th e limiting, T — > 0, uniform 
susceptibility. Consequently the ratio xo /tot — constant 
independent of (U/W), within mean-field, differently 
from the paramagnon or Gaussian results which yield 
nix oc Ln\S\, as we saw previously and consequently to a 
diverging ratio Xo/ m T as \S\ — > |jj 0. 

It is interesting to compare these different predictions 
with measurements at low temperatures of the electronic 
and magnetic Gruneisen parameters of bulk Palladium, 
the prototype of an itinerant 3 — d nearly ferromagnetic 
system. These pressure experiments in which the ratio 
(U /W) is varied due to it's dependence on volume V yield 



I 



9Ln [xo/m T ] 



dLnV 



dLnxo dLnnriT 
dLnV ~ dLnV 



-2.8 - 2.2 



-0.6 ±0.5 



This near cancellation |J implies Xo/ m T ~ constant in 
good agreement with the mean-field scaling prediction. 
The analysis above points out that the true critical be- 
havior of the renormalized, nearly magnetic, Fermi liq- 
uid below the coherence line is governed by mean-field 
exponents. Since the existence of paramagnons in nearly 
ferromagnetic systems is well established |l0) we expect 
that these Gaussian type of fluctuations become impor- 
tant above this line as we shall further discuss. 

The breakdown of hyperscaling in the Fermi liquid 
regime below the coherence line can be more specifically 
discussed within Hertz approach to quantum magnetic 
phase transitions 0. In this case it can be seen how it 
arises. It is due to the dangerous irrelevant quartic inter- 
action uq in the field operators, for d+z > 4 m. The scal- 
ing form of the singular part of the T = free energy can 



be written in this approach as / oc \8\ 



F[u \S\- 



3 



§,§. The dan gerous irrelevant nature of u is a conse- 
quence that F[x — > 0] oc This asymptotic behav- 
ior of the scaling function can be easily obtained if we 
consider the expansion of the free energy in the ordered 
magnetic phase in terms of the order parameter 77, i.e., 
/ oc rj 2 + uaT] 4 + . . . and recall that 77 oc (uq)^ 1 / 2 in Lan- 
dau theory Q. Using this result for F[x — > 0] we get 

/ oc |<^ (<i+2:) / (u |£| ±tFl ) which leads to / oc |<5| 2 and 
consequently to the mean-held exponent a = which 
violates Eq.8 for d + z > 4 since v = 1/2. 

The argument presented above, for T = 0, should hold 
below the coherence line, i.e., for T < T c . For T > T c , the 
interaction uo couples to temperature giving rise to a new 
effective field (uqT) which is now relevant Q. Conse- 
quently the previously discussed mechanism for violation 
of hyperscaling associated with the dangerous irrelevance 
of uo is inoperative for T > T c . Then, for T > T c and in 
particular at 5 — 0, there is no reason to expect violation 
of hyperscaling and Gaussian fluctuations should become 
dominant. This is indeed indicated by the experiments 
in heavy fermion systems as we will discuss. 



V. THE HEAVY FERMION FIXED POINT 

The analysis of experimental data on heavy fermions 
below the coherence line, T < T c , has led to the follow- 
ing empirical relations between the exponents associated 
with the zero temperature fixed point: 2 — a — vz and 
<ph = vz p"l| ] , where </>/j controls the scaling of the uniform 
magnetic field @||]. The former relation clearly violates 
hyperscaling, Eq.8, for d = 3 suggesting that the fixed 
point governing the physics of heavy fermions when ap- 
proached from below the coherence temperature (T < T c ) 
is characterized by classical mean-field exponents at least 
in d = 3. As we pointed out above mean-field exponents, 
contrary to the Gaussian ones, violate hyperscaling above 
the upper critical dimension d c , i.e., in the case of quan- 
tum transitions, for d + z > d c ||. Notice that the rela- 
tions 2 — a = vz and <ph — vz, which lead to 2 — a — <\>% 
imply that the T = uniform susceptibility behaves as 
Xo oc \S\ and the thermal mass, mx oc |<5| as the 
critical point is approached at zero temperature. This be- 
havior of xo shows the importance of ferromagnetic fluc- 
tuations and the reason for the enhanced uniform suscep- 
tibility in heavy fermions. Since the staggered suscepti- 
bility also diverges at the transition, the zero temperature 
critical point which governs the physics of heavy fermions 
is in fact a multicritical point. This led us to consider 
that the heavy fermion fixed point is in the universal- 
ity class of the classical tricritical point with a = 1/2, 
v = 1/2, 7 = 1, (f>T = vz/(ph = 1 and furthermore with 
the dynamic exponent z assuming the value z = 3. This 
yields for the coherence line the equation T c oc |(5| 3 ^ 2 . 
Since vz > 1, this line rises from the T — axis with 
zero derivative reducing the Fermi liquid region close to 
the critical point for (J/W) > (J/W) c . This may be re- 



sponsible for the ubiquity of non- Fermi liquid behavior in 
heavy fermions driven to a non-magnetic state by chemi- 
cal or external pressure (ljjl^l . There are of course other 
possible scenarios as the nearly antiferromagnetic Fermi 
liquid m with z = 2 such that vz — 1. For disordered 
systems one finds z = 4 associated with diffusive modes 
B and v = 1/2 such that vz = 2. 



VI. NON-FERMI LIQUID BEHAVIOR 

Fig.l shows the phase diagram of heavy fermions for 
vz > ip = 1 such that the critical line vanishes linearly 
close to the T = critical point as observed || . It is clear 
from this figure and our interpretation of the crossover 
line that a system just at the critical point, i.e., with 
(J/W) = (J/W) c , \S\ = 0, does not cross the coherence 
line and consequently does not enter the Fermi liquid 
regime. In previous papers Q|| we have obtained the 
properties of such a system under the assumption that 
ijj = vz which is the content of the so-called generalized 
scaling hypothesis |{L4) . The experiments however indi- 
cate that this may not hold in heavy fermions and in 
fact vz > ip = 1 B]- We shall then calculate the ther- 
modynamic properties of the system at the critical point 
for the case vz > ip. For this purpose we introduce a 
set of critical exponents, which we distinguish by a tilde 
and describe the singularities along the critical Neel line 
p5| . It is clear that these tilde exponents are differ- 
ent from those of the zero temperature critical point . In 
the renormalization group language this is a consequence 
that temperature is a relevant field and consequently the 
finite temperature transitions are governed by another, 
T/0, fixed point. 



A. Free energy 

Let us consider the expression for the free energy den- 
sity, 



f<x\j-j c (T)f- a F[t] 



where 



and 



t = 



T/J 



\j-Jc(T)\" 



(12) 



(13) 



(14) 



The scaling function F[t) has the following asymptotic 
behaviors. When t — * 0, F[0] = constant. On the 
other hand for j — > j c (T), he., for t — » 00 we should 
recover the critical behavior of the finite temperature 
antiferromagnetic transition, then F[t — > 00] oc t x with 
x = ^— ^ H. This guarantees that in this limit / oc 

\j-j c (T)\ 2 - & A(T) or / cx \T - T N (j)\ 2 ~ & A(T) where 



4 



the amplitude A(T) = T »* . At the critical point, i.e., 
(J/W) = (J/W) c , we obtain 



focT- 



(15) 



Let us consider some particular cases. 

i) ip = vz. In this case / oc T~^~ . Using hyperscaling, 
i.e., Eq.8, we find a correction to the singular part of the 
free energy density / oc which yields for the specific 
heat C/T oc . In particular for d = z this yields 
C/T oc LnT. On the other hand notice that using the 
empirical mean- field relation 2— a = vz we get / oc T and 
the specific heat is determined by analytic contributions 
to the free energy such that C oc T. 

ii) ip < vz. In this case the specific heat is generally 
given by 

d 2 f 



Hyperscaling yields 



(2-sKt 



-4,) + 4, v (d-z) 



d 2 f 



for tp = 1, d = z = 3 and vz = 3/2 wc get C/T oc T~ . 
For z = 2, such that vz = ip = 1, we obtain C/T oc T 1 / 2 
as found previously by Millis using a different approach 
H . In every case we should consider an additional tem- 
perature independent regular contribution to C/T. 

The mean-field result is obtained using the empirical 
relation 2 — a = vz, we get 



C oc T — 

In the case ip — 1, vz = 3/2 and a « we get 
a stronger than logarithmic divergence for the thermal 
mass, namely, C/T oc T -1 / 3 . For ip = 1, vz = 2, which 
is the disordered case, we get C/T oc T -Q / 2 . For <5 small 
this gives essentially C/T oc LnT. 



B. Order Parameter Susceptibility 

We start with the scaling expression 
Xs ^\j~ Jc (T)\-^F(t) 
which close to the critical Neel line becomes 

For {J/W) = {J/W) c we get 

Xs ocT~^"-)~- (16) 

Taking -0 = 1 and 7=1, appropriate for d+z > d c , yields 
% s oc T a - for z^z = 3/2 and \s °c T for vz = 1. 



C. Uniform Susceptibility 

If we take into account that the uniform susceptibil- 
ity does not diverge at the critical Neel line we get the 
following expression for xo at (J/W) c , 



Xo oc T~ 



(17) 



Now using hyperscaling, which is consistent with Gaus- 
sian exponents and the relation tph = vz [ [TT| , we obtain 



Xo oc T = 



(18) 



which, assuming z = 3, leads to a constant or at most 
a weakly logarithmic divergent susceptibility at (J/W) c . 
In this case also dm/ dh oc Lnh or constant. 

On the other hand taking the dynamic exponent z — 
2, appropriate for antiferromagnetic fluctuations, yields, 
Xo oc T 1 / 2 , at 5 = 0, eventually with an additional, tem- 
perature independent, regular contribution. We also get 
dm/dh oc h 1 / 2 in this case. The above temperature de- 
pendence of xo has in fact been observed Q . 

Notice that using the mean-field empirical relation 
2 — a = vz and vz = <j>h we find that at the critical 
point xo oc T _1 which is the expected mean- field result 
for a Neel transition ocurring at T = 0. This singular 
behavior is quite different from that observed. For the 
differential uniform susceptibility mean-field yields at the 
critical point, \h = dm/dh oc hr x . 

The behavior of the uniform susceptibility at \b\ = is 
relevant as concerns the possibility of superconductivity 
at this point. In fact a xo diverging with temperature at 
the critical point would appear detrimental for supercon- 
ductivity since it signals the presence of unscreened net 
magnetic moments. 



D. Nuclear Relaxation Time 



The nuclear relaxation time is given by |ll|: 



±-=k B Tj2\A q 



,Imx(q,ui N ) 



(19) 



where A q is the form factor and Imx is the imaginary part 
of the wavevector and frequency dependent susceptibility 
X {q, u>) of the system at the nuclear resonance frequency 
Wjv". This expression has contribution from all wavevec- 
tors q. Let us consider the contribution of wavevectors 
near q = Qq the wavector characterizing the ordered anti- 
ferromagnetic phase. Using the dynamic scaling hypoth- 
esis 111 



Imx(.q,«J N )=x.(T)D(qZM') 



(20) 



where £ is the correlation length and z the dynamic expo- 
nent and substituting in the previous expression we get 



5 



— = k B T 



d Xs{T)D{qiM z 



d d q 



which yields 



oc T Xs (T)£ 



z — d 



(21) 



(22) 



Note that the quantity £ 2 d scales as the thermal mass, 
i.e., t\ z ~ d oc m T oc C/T such that ^ cx x s (T)C(T). 
Consequently if C/T oc LnT H, we expect to find using 
the previous results for Xs at least a logarithmic diverging 
linewidth, i.e., Jr oc LnT close to the critical point (7 > 1 

0). 

For completeness we give the temperature dependence 
of the correlation length at the critical point, we get, 
£ oc r p-v/il>-\-{i>—v)/v» m 

We now return to the phase diagram of Fig.l. The 
form of the scaling function Eq.12 together with Eqs. 13 
and 14 allow to identify two additional crossover lines in 
this phase diagram, T c \ and T C 2, whenever tp ^ vz 
Both these lines are governed by the same exponent ip of 
the critical Neel line §,§, i.e., T cl oc \jf and T c2 oc \jf 



The line to the left of (J/W) c , T c \, marks the onset of 
the classical regime where the singularities close to the 
Neel line are described by the tilde exponents. In the 
non-critical side of the phase diagram the temperature 
dependence of the physical quantities in the regime be- 
tween T C 2 and the coherence line (T c < T < T^), are 
obtained from the expressions deduced before for S = 
and taking ip — vz since ip is the relevant crossover ex- 
ponent in this temperature interval . In this region we 
expect hyperscaling to hold and that Gaussian fluctua- 
tions become dominant. The thermal mass for example 
in this temperature region (T c < T < T&) is given by 
C/T oc LnT, for z = 3. 



E. Resistivity at the critical point 

The resistivity at the critical point can be easily ob- 
tained generalizing previous calculations that consider 
the scattering of conduction electrons by different type 
of bosons ||18|— p0| . The resistivity can be expressed, us- 
ing the fluctuation-dissipation theorem, in terms of the 
dynamic susceptibility y(k,oj) associated with the ele- 
mentary excitations [|18| . In three dimensions it is given 

by [pa 



Pc = A(3 



du- 



1] [1 



,-/3w] 



21-, 



dkk 3 \F k \ 2 Imx{k,uj) 



(23) 



where A — 



iehKi' 



(3 = -£jpT an d Fk is a form factor. A 
similar expresion but with the integral in k extending in 
the interval [0, +00] and weighted by a factor k 2 rather 



than fc 3 is obtained if the scattering of the electrons by 
the bosons does not conserve momentum f20[| . This may 
occur in a disordered alloy in which case the boson prop- 
agator should be averaged over all wavectors |2l]] . In this 
case there is always a finite residual resistivity. In order 
to make contact with previous approaches Jl8|,[l9| we con- 
sider, without loss of generality, an explicit expression for 
the imaginary part of the dynamic susceptibility at the 
critical point 



Imx(k,Uj) = Xstatic — 



(24) 



where Q = {u/k l ), A = Dk m . Taking F k = Ck p and 
making a double change of variables in the integrals for 
the resistivity we obtain pc| ] 



p oc T 



and in the case momentum is not conserved 



(25) 



Ap 



Alloy 



oc T~ 



where z — I + m is the dynamic critical exponent. We 
neglect the k-dependence of the form factor (p = 0) and 
the temperature dependence of Xstatic- For p = 0, I = 1 
and z = 3 this gives the result p oc T~s obtained by 
Mathon due to critical paramagnons. In the same 
case but without momentum conservation one finds at 
the critical point Ap AUoy oc T~s . For diffusive modes 
with I = 2 and z = 4 we get, p oc T 3 / 2 and without 
momentum conservation Ap AUoy oc T 5 / 4 at the critical 
point. 

The case of antiferro-paramagnons, I — 0, z — 2 yields, 
p oc T 2 and Ap AUoy oc T2. In the former case, taking 
into account the temperature dependence of Xstatic — 
Xs oc T~ x obtained before, we get a linear temperature 
dependent resistivity at the critical point as observed [|| . 

It is interesting to consider the high temperature be- 
havior (T >> T c ) of Eq. 23 for the resistivity away from 
the critical point. Using the Kramers-Kronig relation we 
find 



P k Ak B T 



2k r 



dkk 3 \F k \ Xstatic(k) 



This high temperature linear behavior is characteristic of 
alloys where a crossover to a T 2 behavior at low temper- 
atures is also observed [E2|. 



VII. DISCUSSION AND CONCLUSIONS 

Using an expansion of the renormalization group equa- 
tions close to the zero temperature fixed point describing 
the physics of heavy fermions we have obtained a critical 
line which vanishes linearly close to the zero temperature 
critical point independently of the crossover exponent vz 



G 



| p3| . We have then generalized our previous scaling ap- 
proach for the case the shift exponent describing the crit- 
ical Neel line is different from the crossover exponent of 
the coherence line. In particular we have obtained the 
properties of the non-Fermi liquid system at the critical 
point, i.e., for (J/W) — (J/W) c . In the general case that 
'tjj =/= vz we have to consider two new crossover lines in the 
phase diagram of Figure 1 , both governed by the shift ex- 
ponent if) |g[. In the region between T c \ and T C 2, which 
includes the non-Fermi liquid regime above the critical 
point, the tilde exponents which describe the singulari- 
ties along the Neel line appear explicitly in the tempera- 
ture dependence of the physical quantities besides those 
associated with the T = fixed point. We have con- 
sidered different scenarios or universality classes for the 
T = transition. 

Our previous analysis of the pressure dependence of 
several physical quantities in different heavy fermion sys- 
tems below the coherence temperature led to the rela- 
tions 2 — a = vz and (f>h — vz. The former equality 
is inconsistent with hyperscaling implying that in this 
region of the phase diagram hyperscaling is violated and 
the physical properties are described by mean-field expo- 
nents. Due to the generality of the concept of coherence 
line this result is also valid for nearly ferromagnetic sys- 
tems as our analysis of the pressure dependence of the 
Wilson ratio of Pd has shown. 

While mean-field exponents are appropriate to de- 
scribe the Fermi liquid regime for T < T c , above the 
coherence line and in particular at the critical point hy- 
perscaling is eventually restored. Experimentally this is 
shown, in the case of heavy fermions, by the saturation of 
the uniform susceptibility at very low temperatures and 
the weak divergence of the thermal mass at the critical 
point. This can be obtained using hyperscaling while ac- 
cording to mean-field %o is expected to diverge strongly 
(Xo c< T _1 ) and C'/T has a stronger than logarithmic di- 
vergence. In the context of the phenomenological scaling 
approach it is not possible to determine a priori the 
validity of hyperscaling in the different regions of the 
phase diagram. However using Hertz approach to quan- 
tum magnetic phase transitions we have discussed how 
a mechanism for the breakdown of hyperscaling at low 
temperatures, i.e., below the coherence line, in fact be- 
comes ineffective for T > T c and hyperscaling is restored 
above this line. 

We have pointed out the relevance of ferromagnetic 
fluctuations in heavy fermions as indicated, for example, 
by the enhanced uniform susceptibility. These fluctua- 
tions, which are associated with the dynamic exponent 
z = 3, may be due to the existence of ferromagnetically 
correlated planes which in turn order antiferromagneti- 
cally. This type of correlations is often found in metam- 
agnetic systems |24| . Our analysis has shown that in the 
case vz ^ ip we do not recover in general, for z = 3, the 
result C /T cx LnT at S = but away from the critical 
point, i.e., for T c < T < T C 2- This region extends down 



to very low temperatures close to 5 = since T c cx |<5| 3 ^ 2 
making it difficult to identify the critical point. Also 
away from 5 — and T > T c , it is easier to obtain a 
linear temperature dependent resistivity as outlined in 
section 6.5. 

We found that taking the dynamic exponent z = 2 
yields the observed temperature dependence of the uni- 
form susceptibility at 8 = ||. In this case, even with- 
out considering the existence of hot lines in the Fermi 
surface , we could obtain a linear temperature depen- 
dent resistivity at 5 = [ p6[ . It is clear from our previous 
remarks that the case of z = 2 is also the most favorable 
for the appearance of superconductivity at the critical 
point {27|. 

We hope further experiments are carried out to dis- 
tinguish between the different scenarios and determine 
unambiguously the universality class of the zero temper- 
ature heavy fermion fixed point. The general scaling the- 
ory formulated above provides a powerful tool to acom- 
plish this task. 

The scaling theory has proved to be useful to describe 
the properties of Ce based heavy fermions. However in 
the case of U compounds, probably due to the exten- 
sion of the 5 — / orbitals, new phenomena appear as the 
coexistence of long range magnetic order and supercon- 
ductivity |2§| ] which would be interesting to incorporate 
in the present approach. On the other hand there is 
a family of materials based on Uranium which exhibit 
many of the non- Fermi liquid properties discussed above 
p9| p0[ . Whether these systems, which have high resid- 
ual resistivities, may be described by the ideas discussed 
here, possibly with differents exponents or within a single 
impurity scenario J3~i[ | remains to be investigated. 

ACKNOWLEDGMENTS 

I would like to thank Profs. F. Steglich, H. v. Lohney- 
sen, A. Loidl, Raimundo R. dos Santos and Drs. E. Mi- 
randa, N. Boncstccl, and V. Dobrasavljevic for useful dis- 
cussions. I am particularly thankful for Dr. E. Miranda 
for a critical reading of the manuscript. I also thank 
CNPq-Brasil for partial financial support. 



[1] Continentino M. A., Japiassu G. M., Troper A.: Phys. 
Rev. B39 , 9734 (1989); Continentino M.A.: J.de 
Physique I 1 , 693 (1991). 

[2] Lohneysen H.v., Pietrus T., Portisch G, Schlager H.G, 
Schroder A., Sieck M., Trappmann T.: Phys. Rev. Lett. 
72 , 3262 (1994); Bogenberger B., Lohneysen H.v.: 
Phys.Rev.Lett. 74 , 1016 (1995); Steglich F., Geibel G, 
Gloos K., Olesch G, Schank G, Wassilew G, Loidl A., 



7 



Krimmel A., Stewart G.R.: J. Low Temp.Phys., 95 , 3 
(1994). 

[3] Contincntino M.A.: Phys.Rev. B47 , 11587 (1993); ibid 

Phys.Rep. Vol.239, 179 (1994). 
[4] Hertz J.: Phys.Rev.B14 , 1165 (1976). For an early 

application of zero temperature scaling ideas to heavy 

fcrmions see: Lawrence J.: J.Appl.Phys.53, 2117 (1982). 
[5] Pfeuty P.: J.Phys.C: Solid State Phys. 9 , 3993 (1976). 
[6] Millis A.J.: Phys.Rev. B48 , 7183 (1993); Zulicke U., 

Millis A.J.: Phys.Rev. B51 , 8996 (1995). 
[7] Doniach S., Engelsberg S.: Phys.Rev.Lett. 17, 750 

(1966). 

[8] Toulouse G., Pfeuty P.: Introduction au Groupe de 
Renormalization et a ses Applications, Grenoble, Presses 
Universitaires de Grenoble, 1975; Fisher M.E. in Crit- 
ical Phenomena, ed. F.J.W.Hahne; Proa, Stellenbosch, 
South Africa (Springer, Berlin 1983) p.l. 

[9] Gerhardt W., Razavi F., Schilling J.S., Huser D., My- 
dosh J.A.: Phys.Rev.B24 , 6744 (1981); Barron T.H.K., 
Collins J.G., White G.K.: Advances in Physics 29 , 609 
(1980). 

[10] Levin K., Vails O.T.: Phys. Rep. Vol.98, 1 (1983). 

[11] For an isotropic antiferromagnet at 8 — this has been 

obtained theoretically by Sachdev S., Z.Phys. B94 , 469 

(1994). 

[12] Budko S. et al.: to be published in Physica B . 
[13] Loidl A. et al.: private communication. 
[14] Pfeuty P., Jasnow D., Fisher M.E. Rhys. Rev. BIO , 2088 
(1974). 

[15] The tilde exponents are most probably those of the 3 — d 

Heisenberg model. 
[16] Bourbonnais C, Wzieteek P., Creuzet F., Jerome D., 

Batail P., Bechgaard K.: Phys.Rev.Lett. 62 , 1532 

(1989). 

[17] Halperin B.I., Hohenberg P.: Rev.Mod.Phys. 49 , 435 
(1977). 

[18] Mills D.L., Lederer P., J.Phys.Chem.Sol. 27 , 1805 
(1966). 

[19] Mathon J., Proc.Roy.Soc. A306 , 355 (1968). 

[20] Rivicr N., Mensah A.E., Physica B91 , 85 (1977). 

[21] Continentino M.A., Rivier N., J.Phys.F: Metal Physics, 

8 , 1187 (1978). 
[22] Kaiser A.B., Doniach S., Intern. J. Magnetism, 1 , 11 

(1970). 

[23] Notice that if the T — transition is first order then 
dT N /dP -c -oo as T N -> or P -> P c where 
P c is the critical pressure; Samara G.A., Terauchi H., 
Phys.Rev.Lett. 59, 347 (1987). 

[24] Haen P., Besnus M.J., Mallmann F., Mignot J.-M., Lejay 
P., Meyer A.: Physica B 199 & 200 , 522 (1994). 

[25] Hlubina R., Rice T.M., Phys.Rev. B51 , 9254 (1995). 

[26] This has also been obtained by Aliev F.G., JETP Lett. 
61, 1017 (1995), using a different approach. 

[27] Movshovich R., Graf T., Mandrus D., Thompson J.D., 
Smith J.L and Z.Fisk, to be published; Grosche F.M., 
Julian S.R., Mathur N.D. and Lonzarich G.G to be pub- 
lished. 

[28] Steglich F., Geibel C, Modler R., Lang M., Hellmann 
P., Gegenwart P.:J.Low Temp. Phys. 99 , 267 (1995); 
Steglich F., Buschinger B., Gegenwart P., Geibel C, Hel- 



frich R., Hellmann P., Lang M., Link A., Modler R., Jac- 
card D., Link P. Proceedings of the International Confer- 
ence on Physical Phenomena in High Magnetic Fields II, 
Tallahassee (1995) to be published. 

[29] Andraka B., Physica B 199 & 200 , 239 (1994). 

[30] Tsvelik A.M., Reiser M., Phys.Rev.B48 , 9887 (1993). 

[31] Dobrasavljevic V., Kirkpatrick T.R., Kotliar G., Phys. 
Rev.Lett. 69, 1113 (1992). 



8 



